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^ ; ABSTRACT 
> 

Observations by the RHESSI satellite of large polarization of the prompt 7-ray 
emission from the Gamma Ray Burst GRB021206 (Coburn & Boggs 2003) imply that 
■ the magnetic field coherence scale is larger than the size of the visible emitting region, 

~ R/T, where R is the radius of the flow, T is the associated Lorentz factor. Such 
fields cannot be generated in a causally disconnected, hydrodynamically dominated 
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outflow. Electromagnetic models of GRBs (Lyutikov & Blandford 2002), in which large 
scale, dynamically dominant, magnetic fields are present in the outflow from the very 



O 

beginning, provide a natural explanation of this large reported linear polarization. We 
derive Stokes parameters of synchrotron emission of a relativistically moving plasma 
with a given magnetic field configuration and calculate the pulse averaged polarization 
fraction of the emission from a relativistically expanding shell carrying global toroidal 
magnetic field. For viewing angles larger than 1/T the observed patch of the emitting 
shell has almost homogeneous magnetic field, producing a large fractional polarization 
(56% for a power-law energy distribution of relativistic particles dn/de oc e~ 3 ). The 
maximum polarization is smaller than the theoretical upper limit for a stationary plasma 
in uniform magnetic field due to relativistic kinematic effects. 

Subject headings: gamma rays: bursts — MHD — polarization 



1. Introduction 

Origin of magnetic fields in Gamma Ray Bursts (GRBs) is one of the central unresolved issues. 
In the standard fireball scenario (e.g. Piran 1999, Meszaros 2002 and references therein), magnetic 
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field does not play any dynamical role. The near-equipartition field invoked in the emission region 
is assumed to be generated locally at relativistic shocks by plasma instabilities (e.g. Medvedev & 
Loeb 1999). Initially, the spatial scale of such fields is microscopically small, of the order of the ion 
skip depth, 6 ~ c/uj Pi i (co p j is the ion plasma frequency). Though the typical scale of magnetic field 
fluctuations may grow due to inverse cascade, even in the unlikely case that such growth proceeds 
at the speed of light the resulting polarization is expected to be smaller than 10% (e.g. Gruzinov 
& Waxman 1999). 

The recent detection by the RHESSI satellite of large polarization in the prompt 7-ray emission 
(Coburn & Boggs 2003) places severe constraints on the GRB models. It implies that magnetic field 
coherence scale is larger than the size of the visible emitting region, ~ R/T, where R is the distance 
form the center and T is a bulk Lorentz factor of relativistically expanding emission region. Such 
fields cannot be generated in a hydrodynamically dominated outflow, which is causally disconnected 
on large scales. Thus, the large scale magnetic fields should be present in the outflow from the very 
beginning. In fact, as we argue below, such fields must be dynamically dominant, carrying most of 
the energy of the outflow. 

Building upon earlier models of electromagnetic explosions (e.g. Usov 1992; Thompson 1994; 
Smolsky & Usov 1996; Meszaros & Rees 1997), Lyutikov & Blandford (2002, 2003) developed an 
electromagnetic model of GRBs which assumes that rotating, relativistic, stellar mass progenitor 
(e.g. "millisecond magnetar", Usov 1992) loses much of its spin energy in the form of an electro- 
magnetically dominated outflow. A stellar mass relativistic progenitor is born with angular velocity 
O ~ 10 4 s _1 and dynamo generated magnetic field of B s ~ 3 x 10 14 G. Then the total rotational 
energy, E ~ IQ? /2 ~ 5 x 10 52 erg (for a 1.4 M object) is available to power GRB bursts, while 
the dipole spin-down luminosity Lp ~ i? 2 r^$7 4 /c 3 ~ 10 49 ergss _1 is about the luminosities of cos- 
mological 7-ray bursters. In this model the energy to power the GRBs comes eventually from 
the rotational energy of the progenitor. It is first converted into magnetic energy by the dynamo 
action of the unipolar inductor, propagated in the form of Poynting flux dominated flow and then 
dissipated at large distances from the sources. 

A rapidly spinning magnetar with a complicated field structure will form a relativistic outflow. 
We suggest that magnetic field in the wind quickly rearrange to become predominantly axisym- 
metric. There is a good precedent for this behavior in Ulysses observations of the quiet solar wind 
(McComas et al. 2000) which reveal that, despite the complexity of the measured surface magnetic 
field, the field in the solar wind quickly rearranges to form a good approximation to a Parker (1960) 
spiral. The situation in the far field will then resemble that first analyzed by Goldreich & Julian 
(1969) and the characteristic scale length in the far field is the cylindrical radius from a polar axis, 
rather than the wavelength. 

During the relativistic expansion, most of the magnetic energy carried by axisymmetric toroidal 
magnetic fields is concentrated in a thin shell with thickness Ar ~ ct s ~ 3 x 10 12 cm inside 
a contact discontinuity separating the ejecta from the shocked circumstellar material. At the 
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contact discontinuity toroidal magnetic field balances the ram pressure of the circumstellar material 
Bff, ~ 4r 2 \j Tip cy ±(? ■ Both B^ and V depend on the angle between a given point on the shell and 
the polar axis, defined as the axis of rotation of the progenitor. This results in a non-spherical, 
relativistic expansion of the shell. In particular, for laterally balanced expansion T ~ 1/ sin#. The 
current carrying shell becomes unstable due to development of the current driven instabilities at a 
radius ~ 10 16 cm. This leads to acceleration of pairs which emit 7-rays by synchrotron radiation. 

A distinctive feature of the electromagnetic model is that the causal connection is better than 
in the hydrodynamic models. Initially, close to the central source the subsonic flow is fully causally 
connected. As the flow is accelerated by magnetic (and partially by pressure) forces it becomes 
supersonic, strongly relativistic and causally disconnected over small polar angles A6> ~ 1/T. Later, 
magnetically dominated flows quickly reestablish causal contact over large polar angles and become 
fully causally- connected again after a time t c ~ t s r 2 , where t s ~ 100 sec is the source activity life 
time. This is drastically different from hydrodynamic flows which remain causally disconnected over 
polar angles larger than 1/I\ Thus, during expansion the causal behavior of the flow resembles the 
behavior of cosmic fluctuations during inflation: as the flow expands, angular scales ~ T _1 "enter 
the horizon" re-establishing causal contact that was lost during acceleration. 

To illustrate this behavior, consider propagation of a sound-type disturbance emitted by a 
point source located on the relativistically moving shell at radius R cm - Let the typical signal speed 
in the plasma rest frame be (3 S . In appendix A we show that for sub-Alfvenic ejecta (magnetically 
dominated flows can be strongly relativistic, but still sub-Alfvenic!) a relativistically expanding shell 
re-establishes a causal contact over the visible patch of 1/T in just one dynamical time scale (after 
doubling in radius). If the ratio of the magnetic to particle energy density in the cold, magnetized 
plasma is a = ub/u p ^> 1 (Kennel & Coroniti 1984), then the Alfven (and fast magneto-sound) 
velocity is Pa = C\J a/ (1 + a). The requirement that the expansion velocity be sub-Alfvenic then 
implies that cT < u\ or a S> 1. Therefore, the condition that magnetic fields have a coherence 
scale larger than R/T requires that the magnetic fields be energetically dominant in the flow. 

Hydrodynamic (e.g. fireballs, Piran 1999, or external shocks, Dermer 2002, with a <C 1) or 
hydromagnetic models (a ~ 1, e.g. Spruit et al. 2000, Drenkhahn & Spruit 2002, Vlahakis & Konigl 
2003) could also have a large scale ordered magnetic field (cannonballs, e.g. Dar 2003, need to rely 
on energetically inefficient Compton scattering and strong flow inhomogeneities on the angular 
scale ~ 1/r, Coburn & Boggs (2003), section 3; we consider this prohibitively tight constraints). 
In all cases, the whole outflow is in causal contact close to the source and may have a large scale 
magnetic field which will be carried with the flow. In hydro-dominated models, after the causal 
contact is lost, different parts of the flow cannot communicate and thus will evolve differently, 
depending on the local conditions. Only under strict homogeneity of the surrounding medium 
and of the ejecta the two causally disconnected parts of the flow will have similar properties. On 
the other hand, since magnetically dominated outflows can quickly communicate information (e.g. 
magnetic pressure) over large polar angles, they can have quasi-homogeneous properties despite 
possible inhomogeneities in the circumstellar medium and in ejecta. 
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Assumption of electromagnetically dominated flow must eventually break down, since the 
ejecta need to dissipate magnetic energy to produce high energy emission. In the emission region 
the plasma is expected to be close to equipartition as magnetic field is dissipated to accelerate 
electrons (generally, equipartition is needed for effective emission). But unlike the hydrodynamic 
models where equipartition is reached by amplifying weak magnetic fields at the shocks, in the elec- 
tromagnetically dominated model the equipartition is reached by dissipation of initially dominant 
magnetic field, as, for example, happens in solar flares. 

In this paper we calculate the Stokes parameters for the prompt GRB emission emerging in 
electromagnetic model as a function of the viewing angle (angle between the line of sight and the 
polar axis of the flow). We assume that magnetic field in the emission region is dominated by 
the toroidal field and is concentrated in a thin shell AR near the surface of the shell expanding 
with the Lorentz factor T(9). Synchrotron emission is produced by an isotropic population of 
relativistic electrons with the power law distribution in energy. In the present work we calculate 
Stokes parameters averaged over the duration of the GRB pulse, deferring the time-dependent 
calculations to a later work. One expects that the polarization fraction will be maximal in the 
beginning of the pulse, slightly decreasing towards the end as larger emitting volumes become 
visible. We set the speed of light to unity, c = 1, in all the expressions to follow. 



2. Calculation of Stokes parameters 

Consider a quasi-spherical thin emitting shell (Fig. 1) viewed by an observer. Below we denote 
all quantities measured in the local frame comoving with an emitting elementary volume with 
primes, while unprimed notations refer to the quantities measured in the explosion frame. Let 
r, 9, and 4> be the spherical coordinates in the coordinate system centered at the center of the 
shell, and x, y, and z be the rectangular coordinates with the origin at the center of the shell. 
The symmetry axis of the shell is the z-axis. The toroidal magnetic field in the shell is in the 
0-direction. The observer is located in the x—z plane. The components of all vectors written below 
are the components with respect to the rectangular coordinate system x, y, and z. The shell 
expands quasi-spherically with an angle-dependent Lorentz factor T(9). An element of the shell 
moving radially with the velocity v = /3{sin 9 cos 4>, sin 9 sin 4>, cos 9} emits a burst of synchrotron 
radiation in the direction of unit vector n = {sin# b, 0, cos 9 Q b} when viewed in the observer frame. 

Several key ingredients need to be taken into account (e.g. Cocke & Holm 1972, Blandford & 
Konigl 1979, Bjornsson 1982, Ginzburg 1989). First, the synchrotron emissivity depends on the 
direction between the emitted photon and the magnetic field in the plasma rest frame. Second, as 
the emission is boosted by relativistic motion of the shell, the position angle of the linear polarization 
rotates in the n — v plane. 1 The fractional polarization emitted by each element remains the 
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same, but the direction of polarization vector of the radiation emitted by different elements within 
a visible shell is rotated by different amounts. This leads to effective depolarization of the total 
emission. The theoretical maximum polarization fraction for homogeneous field can be achieved 
only for uniform plane-parallel velocity field. Third, integration along the line of sight (and over 
the emitting solid angle for unresolved sources) is better carried out in the laboratory frame, in 
order to take correct account of the arriving times of photons. 

We assume that the distribution function of emitting particles in the frame comoving with an 
element of the shell is isotropic in momentum and is a power law in energy 

dn = K c e- p dedVdn p . (1) 

Here dn is the number of particles in the energy interval e, e + de, dV is the elementary volume, 
dfl p is the elementary solid angle in the direction of the particle momentum p, K c = K e (r), 
p = constant. 

In this paper we are interested in the polarization structure of the time integrated pulse of 
the emission, and not in its temporal properties. Hence, for simplicity we approximate the 7-ray 
emissivity of the shell as a flash at some time to i n the explosion frame, lasting for At <C AR/c, 
where AR <C Ro is the thickness of the shell at the moment to and Rq is the radius of the shell at 
the moment to- More complicated emission profiles may be easily accommodated. In addition, we 
integrate over the observer time to get an average polarization of the pulse deferring time dependent 
calculations to a later work. 

We also assume that the emission is optically thin and neglect possible plasma propagation 
effects (e.g., depolarization of radiation due to internal Faraday rotation by low energy electrons). 
Since the emitting particles are ultra-relativistic and we neglect conversion of linear to circular 
polarization in plasma, we do not have circular polarization in our model (Stokes V = 0). We 
also neglect a possible tangled component of the magnetic field present in the emission region. We 
assume that the emission originates in a geometrically thin layer AR <C R with the thickness AR 
independent on 9 and neglect variation of the magnetic field and velocity across the layer. Given 
these assumptions, our estimates provide an upper limit on the possible polarization. 

Time-integrated Stokes parameters are calculated in appendix B (Eqns. (B7)). Due to cylin- 
drical symmetry of the model the Stokes parameter U integrates to zero, so that the observed 
averaged polarization fraction is 

_\Q\ _ p+1 fsm6ded(pV 2+ (P- i y 2 \B'smx'\ ip+1)/2 cos2x 
11 = T = p + 7/3 J sin 9d9d4> p2+( P -i)/2 J B , sin x <\ (p+i)/ 2 • W 

Here B' is the magnitude of the magnetic field B' in the frame of an element of the shell, x' is the 
angle that the line of sight in the frame of an element of the shell, n', makes with the magnetic field 
B', and x is the position angle of the electric field vector in the observer plane of the sky measured 
from some reference direction. Doppler boosting factor is V = 1/T(1 — n • v). For toroidal magnetic 
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field Q > 0, so that the observed polarization vector is always along the projection of the flow axis 
on the plane of the sky. 

Evaluation of different quantities in Eqn. (2) is an involved exercise in Lorentz transformations. 
We assume that, in the shell frame, the magnetic field is purely toroidal 

B' = &0(r,0)B' = &4- sin 0, cos 0,0}, (3) 

where B' is the unit vector along B' in the radiating element frame and is the magnitude of the 
field. A photon propagating along the unit vector n in the explosion frame is emitted along the 
direction with the unit vector n' in the radiating element frame: 



r(l-(n-v)) ' u 

Note, that n', n and v lie in the same plane. The angle x' between the photon and magnetic field 
in the radiating element frame is 

(B'.n)+r(B'.v)( r £ I (n.v)-l) 
«*x' = B'.n< = r( l-(n v)) 1 < 5 > 



which gives 

r 2 (i-/^) 2 



sm x = 1 - T^Ti I 6 ) 



where (i = cos 6 cos 6> b + sin 9 sin # b cos </>. 

We also need to evaluate angle x. between a given direction in the observer plane and the 
polarization vector. This is not trivial since polarization vector emitted by each element will 
experience rotation during Lorentz transformation from the shell frame to the laboratory frame 
(Cocke & Holm 1972; Blandford & Konigl 1979). Rotation of the polarization vector is due to 
the rotation of the wave vector in the plane containing vectors n, n', and v, and the requirement 
that the electric field of the wave remains orthogonal to the wave vector. Since wave vectors of 
emitted waves experience rotation by angles of the order of unity, this effect would lead to effective 
depolarization of emission from a medium with non-uniform velocity field even from a homogeneous 
magnetic field. In appendix C we derive general relations for the Lorentz transformation of the 
polarization vector. 

We choose to measure angle x clockwise from the direction parallel to the projection of the 
axis of the flow on the plane of the sky. The unit vector in this direction is {— cos o b, 0, sin # G b}- 
We find (appendix C) 

(1 — 8 a) cos 6 — 3 sin 8 sin o b sin 2 6 

cosx~ 



V(l - (i^) 2 ~ sm 2 #ob sin 2 <f>/T 2 



. _ sin cos 6* — cos# b) 

Sill V \i) 

v / (l-/3/i) 2 -sin 2 ob sin 2 JjT*' 
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In the ultra-relativistic limit, r > 1, general relations simplify and it becomes possible to 
determine analytically the maximum polarization fraction for a given velocity field in the limit 
o bT 3> 1 (appendix D). For p = 3 we find II = 9/16, in excellent agreement with numerical 
calculations. 

On Fig. 2 we plot the map of the polarized emissivity from the shell moving with a constant T 
and with constant as it is seen by the observer on the plane of the sky. I and s are rectangular 
coordinates on the plane of the sky centered at the projection of the center of the shell. Axis s is 
directed parallel to the projection of the axis of the shell on the sky. I and s are normalized such 
that the projection of the shell radius, R s , is a circle of radius 1 in the l-s plane. Thus, 

I = -^-1 • r = sin 6* cos 0, s = — — r • (1 x n) = — cos 6* b sin 8 cos <f> + sin6* b cos 6, (8) 

R s Us 

where 1 = {0,1,0} is the unit vector along I. The arrows on the plots in Fig. 2 are directed 
perpendicular to the unit vector in the direction of the electric field of the wave, e, so that in the 
non-relativistic limit, T 1, the arrows are aligned with the magnetic field B. The length of the 
arrows is proportional to the synchrotron emissivity from the unit volume, i.e., to the expression 
under the integral for I in Eq. (B7). Actual observed intensity is modified by the geometric factor 
proportional to the path of the ray inside the volume of the shell. For R s — r^> AR, the geometric 
factor is 1/fi. For r> 1, Doppler boosting leads to the small effective emitting area of the shell: 
I < 1/r and s < 1/r. Relativistic swing of polarization vector is also clearly visible in Fig. 2. Each 
patch of the shell emits radiation with the same polarization degree, II max = (p + l)/(p + 7/3). 
Due to summation over the areas of the shell with different directions of x, resulting polarization 
degree becomes smaller than II max . 

We are now in a position to estimate polarization fraction (2) integrating Stokes parameters 
(B7) over an expanding relativistic shell. Results for II are shown in Fig. 3. The parameter Q is 
zero for # b = and is small for # Q b < 1/r, because the polar axis falls within the visible patch 
in this case. The magnetic field changes its direction within the visible patch and the resulting 
polarization is reduced. The degree of polarization reaches a limiting value of tens percent when 
observation angle is larger than 1/r. 



3. Discussion 

Large scale ordered magnetic fields produced at the central source provide a simple explanation 
of the recent observations of highly polarized GRB prompt emission by the RHESSI satellite 
(Cob urn & Boggs 2003). In order to retain the coherence of the magnetic field on scales larger than 
the visible patch, ~ R/T, the ejecta must be electromagnetically dominated. The electromagnetic 
model suggested by Lyutikov & Blandford (2002, 2003) provides a solution to the puzzle of how to 
produce large coherent magnetic fields and how to launch a blast wave that extends over an angular 
scale >> r _1 and where the individual parts are out of causal contact. In the electromagnetic 
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model, the magnetic fields are present in the outflow from the very beginning, and the energy is 
transferred to the blast wave by a magnetic shell which is causally connected at the end of the 
coasting phase. 

To prove this point we first found general relations for transformation of polarization direc- 
tion of synchrotron emission produced by relativistically moving source with a given magnetic field 
structure and calculated Stokes parameters for the time averaged synchrotron emission for a par- 
ticular case of relativistically expanding shell containing toroidal magnetic fields. We find that for 
observing angles satisfying o b > 1/T a large polarization fraction IT ~ 60% may be observed (the 
actual spectrum was not measured for GRB021206). The position angle of the polarization is fixed 
by the projection of the progenitor axis on the plane of the sky and thus should not change during 
the burst. 

Another potential source of polarization could be Compton scattering of unpolarized 7-rays. 
If unpolarized 7-rays are initially beamed into a small-angle jet and are scattered by surrounding 
gas, then polarized scattered 7-rays would be distributed nearly isotropically. This will require a 
much higher energy in the initial 7-ray jet than the energy necessary if the narrowly beamed jet is 
observed directly (Coburn & Boggs 2003), thus putting a much tougher requirement for the total 
energy budget of GRB. Compton scattering by relativistically moving wide angle envelope could 
also occur. In this case kinematics of the scattering is similar to the kinematics of the synchrotron 
emission considered in the present work. Therefore, the energetic requirements are also similar to 
the synchrotron mechanism. However, Compton scattered photons do not have preferred direction 
of polarization, which is set by the large scale magnetic field in the synchrotron case. The net 
polarization of scattered photons from uniform spherical shell averages to zero. High polarization 
can be observed only if the shell parameters (T or electron density) vary significantly on the angular 
scale ~ 1/r. We are coming back to highly collimated flow. Therefore, we conclude that Compton 
scattering cannot account for the high degree of polarization of 7-rays emerging from a wide angle 
expanding flow. We note, that synchrotron mechanism results in the electric vector of polarized 
emission directed parallel to the axis of the flow, while the scattered 7-rays would be polarized in 
the direction perpendicular to the axis of flow. 

Several natural correlations between GRB polarization and other parameters follow from the 
model and can be tested with future observations. First, polarization fraction should decrease from 
the beginning to end of the pulse as larger areas of the emitting shell become visible to the observer. 
Second, the maximum amount of polarization is related to the spectrum of emitting particles, being 
higher for softer spectra. This points to a possible correlation between the amount of polarization 
and hardness of the spectrum. 

Our treatment of the prompt emission may also be related to polarization of afterglows. If 
the field from the magnetic shell may be mixed in with the shocked circumstellar material (similar 
to the so called flux transfer events at the day side of Earth magnetosphere), then a comparably 
large fractional polarization may be observed in afterglows as well. In addition, since the preferred 
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direction of polarization is always aligned with the flow axis, the position angle should not be 
changing through the afterglow (if polarization is observed both in prompt and afterglow emission 
the position angle should be the same). Also, polarization should not be related to the "jet break" 
moment. This is in a stark contrast with the jet model, in which polarization is seen only near the 
"jet break" times and the position angle is predicted to experience a flip during the "jet break" 
(Sari 1999). For the same reason large average polarization cannot be due to particular viewing 
geometry, as suggested by Waxman (2003). Current polarization data show that in virtually all 
cases position angle remains constant (Covino et al. 2003a, b; Barth et al. 2003; Bersier et al. 2003, 
see though Rol et al. 2003), while the amount of polarization does not show any correlation with 
the "jet break". This is consistent with the presence of large scale ordered magnetic fields in the 
afterglows. (A model of Rossi et al. (2002) of structured jets also predicts constant position angle, 
but since no large scale magnetic field is assumed the polarization features are still related to the 
jet break times). 

In our calculations we have neglected a random component of the magnetic field which must be 
present in the emission region. Lyutikov & Blandford (2002, 2003) suggested that 7-ray emitting 
electrons are accelerated by current instabilities, somewhat similar to solar flares. Development of 
current instabilities should be accompanied by dissipation of magnetic fields and destruction of the 
magnetic flux. These will generally add a random component to the ordered magnetic field and will 
lead to a decrease in polarization (Korchakov & Syrovatskii 1962). The corresponding calculations 
are in progress. 

An alternative model of GRBs that can feasibly give large scale magnetic fields in the prompt 
emission region is the plerion model (Konigl & Granot 2002; Inoue et al. 2003, see also Lyutikov 
2002), which initially was suggested for afterglows, but may also be extended to include the prompt 
emission by external shock wave (Dermer 2002). In this case, the large scale equipartition magnetic 
fields are created ahead of the expanding GRB ejecta by the preceding explosion of the "supranova" 
(Vietri & Stella 1999). Still this type of models faces similar causality/efficiency problem: if the 
plerion plasma is only at equipartition, a ~ 1, it may be expected to be inhomogeneous on R/T 
scale; if it is strongly magnetized, a S> 1, then the shocks will be only weakly dissipative. 

Implications of the RHESSI results, that GRB flows are electromagnetically-driven, may pro- 
vide an important clue to the dynamics of other astrophysical sources like pulsars, (micro) quasars 
and AGNs. It is quite plausible that all these sources produce ultra-relativistic magnetically dom- 
inated outflows with low baryon density (Blandford 2002). The flow evolution in all these systems 
may proceed in a similar way. Energy, transported primarily by magnetic fields, is dissipated far 
away from the source due to development of current instabilities. Particles are accelerated in local- 
ized current sheets by DC electric fields and/or electromagnetic turbulence producing bright knots 
(in AGNs) and a variety of bright spots in pulsar jet, best observed in the Crab. 
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After the submission of the paper additional observational details of the burst GRB021206 
became available. The spectral index of the burst is a ~ 0.6; corresponding particle index is 
p = 2.2 (Hajdas, private communication). The maximum polarization in our model is then ~ 45%. 
In addition, there were several theoretical developments. Higher polarization (up to 100%) can be 
achieved if the particle distribution is not isotropic (Lazatti, private communication). Granot (2003) 
have performed polarization calculations similar to ours, taking into account a random component 
of the magnetic field. He reached a similar conclusion that it is substantially easier to produce the 
polarization observed in GRB 021206 from an ordered magnetic field. A small discrepancy between 
the results is explained by the difference in the duration of the emission: while we assumed that 
the shell emits during time At <C Ar, Granot (2003) assumed that the emission is more prolonged, 
At > ArT 2 . 
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A. Causal structure of relativistic magnetized outflows 

In this appendix we consider the casual structure of the relativistically expanding magnetized 
shell. We wish to answer the question: "If a surface of relativistically expanding magnetized shell 
is perturbed at a given radius R em and zero polar angle, which points on the surface of the shell 
will be affected after time i" . 

Consider propagation of a sound-type disturbance emitted by a point source located on the 
surface of relativistically moving shell at radius Rem- Let the signal speed in the plasma rest 
frame be [3 S . For simplicity we assume that the shell is moving with constant velocity. If in the 
plasma frame a wave is emitted in the direction 9 em with respect to the flow velocity, 2 then in the 
laboratory frame the components of the wave velocity along and normal to the bulk velocity are 

(3 + f3 s COS gem a _ & sin ^gm ( ... 

Ps,lab,\\ ~ l + m COS 9 em ' P "' lab ' ± ~ r(l + PP S COS em ) ' lAij 

The condition that a waves catches the surface of the shell becomes 

t(3 s sin 9 em 



(R cm + t(5) sin 9 



T(l + pp s cos 9 em y 



(R cm + 0) cos 9 = R cm + f + P ; C ° S9 : m t- (A2) 

1 + P/3 S cos 9 cm 



Eliminating 9 em we find 



In case of isotropic relativistic fluid with internal sound speed (3 S = I/a/3, Eq- (A3) implies that 
the maximum angle that sound waves can reach as t — > oo is 9 ~ \J \7H — 2/7 ~ O.67/7. Thus, 
hydro-dominated relativistic plasma remains causally disconnected on scales 9 ~ 1/r at all times. 

In magnetically dominated medium the situation is drastically different. Consider, for simplic- 
ity, cold magnetically dominated plasma. If the ratio of the magnetic to particle energy density in 
the flow istT = ns/u p »l (Kennel & Coroniti 1984), the Alfven velocity is (3a = a/c/(1 + a). 
Eq. (A3) then becomes 



This implies that two points on the surface of the shell separated by an angle 9 ^ 1 come into 
causal contact after time 

" 'W,/I±*. (A5) 



Rc 



Propagation of waves in a non-uniformly moving medium will generally lead to a change of the wave direction; 
for qualitative estimates we neglect here this effect. This is well justified in the strongly magnetized limit a — > 00 
and/or for small angles 6 < 1. 
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Thus, in strongly magnetized medium, a S> 1, the visible patch of the shell with 9 ~ 1/F re- 
establishes causal contact in one dynamical time t ~ R era /c. 

We may also invert Eq. (A4) to find time needed to establish casual contact over angle 6: 

2i? cm /3rsin| 



t = 



(3 (a - 4(/3 T) 2 sin 2 f (1 + (/3T) 2 sin 2 f )) ' 



x I f(l + «7)(l + (/?r) 2 sin 2 +/3rsin^ ( f3T + 2(1 + (/?r) 2 sin 2 ?-))). (A6) 



For <T < 4T (1 + T ) the maximum causally connected region (for t — > oo) is finite 



. 2 #oo VI + °~ ~ 1 



si »f = ^T)^ (AT) 

For subsonic flow, ^fo > T, 9^ becomes larger than 1/T. For larger a whole shell comes into a 
causal contact after time 

t = flg-^l + I*)) ( 2(/?r)3 + ^Ml + W) + /?r(2 + a)) . (A8) 

It is also instructive to find the emission angle 9 em as a function of time t when an emitted 
wave catches with the surface of the shell: 



C< .S - - - — - nT } , ^ ((i?em + /ft) " V^cm + /?t) 2 1 t/?(2i? cm + /3t)/3 2 ) . (A9) 



PP s (2R em + (3t) 



In an ultra-magnetized plasma (force- free plasma, a — ► oo) relations simplify considerably. 
The region casually connected after time t becomes 



9 11 t 2 
sin - = 



2 -Rem(i?em + ' 

t = 2i? cm r sin r sin °- + ^1 + (/? T) 2 sin 2 |j . (AlO) 

So that points separated by 1/T come into causal contact after t = (1 + V5)-Rcm/2 and the whole 
shell becomes causally connected after i ~ 2i? em r 2 . Emission angle (A9) in the force- free case then 
becomes 

COS ^ m = -2lufw (A11) 
Note that in the plasma rest frame the waves which propagate furthermost in polar angles are 
emitted "backward" in the explosion frame. It is due to this fact why they can "beat" the commonly 
called result that lateral velocity in the laboratory frame cannot be larger than c/T. The latter 
is true only for waves propagating along the surface of the shell (normal to the flow in the shell 
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frame). We can understand then why hydrodynamic sound waves cannot reach large polar angles: 
when emitted "backward", they are advected with the supersonically moving flow. On the other 
hand, in subsonic strongly magnetized plasma fast magnetosound waves can outrun the flow and 
reach large polar angles in the laboratory frame. 

Thus, if the effective signal velocity ( Alfven velocity) in the bulk of the flow is larger than the 
expansion velocity, strongly relativistic magnetized outflows quickly become casually re-connected 
over the visible patch Ad ~ 1/T. For strongly subsonic flows, a » T 2 , points on the shell separated 
by 1/r come into causal contact on a dynamical time scale t ~ R era (in a relativistic hydro dynamical 
flow this never happens). Since this time is fairly short, the global dynamics of the flow is not very 
important which vindicates our assumption of constant expansion velocity. 



B. Pulse-integrated Stokes parameters 

The Stokes parameters are components of the polarization tensor J; s = \ 

Here xi are coordinates in the plane perpendicular to n and there is no circular polarization. The 
pulse-integrated intensity 

=^j dT j dv ^ ( n ' v > r > t + r -^-) • < Bi ) 

where ji s is emissivity, T = t — r cos 0/c is the observer time and integration is over the whole 
emitting region in the explosion frame. 

We approximate emissivity as an instant flash at the moment t = to with the duration At, 
At <C AR/c. We also assume that the whole shell emits uniformly during the flash. Then, the 
emissivity can be expressed as 

j ls {n, u, r, t) = j ls {n, u, t )6(t - to) At [H{r - R ) - H(r - R - AR)} , (B2) 

where 8(x) is Dirac delta function and H(x) is a step function, H(x) = 1 if x > 0, H(x) = if 
x < 0. We first integrate in T keeping all other independent variables (r, 9, (j)) fixed: 

4(f) = jy2 i d(p sinOde / r 2 dr j ts (n, u, r, t ) x 
At [H(r - R ) - H(r - R - AR)} . (B3) 
Taking into account that AR <^ R and integrating in Eq. (B3) over dr, we obtain 

Ji s {u) = AtAR^J^ d<f>J^ sm6d0j ls (n,u,Ro,to). (B4) 

Lorentz transformation of the emissivity to the comoving frame with the element of the shell is 

jls(n,u,t ) =P 2 (n',t )^ s ,(n',P-V) (B5) 
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so we obtain 

?2 i-2-K 



J ls {v) = AtAR^ I"" dcf) ^ smedeV^P-^^j^^n',^). (B6) 
D Jo Jo 

Using synchrotron expressions for j'[, s ,(n', u) in the comoving frame (e.g., Ginzburg 1989) we obtain 



2 t-l-K 



P+l 

?2 /-27T 



Q = n(u)ARAt D 2 °_ 1)/2 ^ # ^ sin^^ 2 +^ 1 )/ 2 | J B'sin X '|(f+ 1 )/ 2 cos2x, 

U = K(v)ARAt — - z) ° 2+{p _ 1)/2 j Q sin^^ 2 +^- 1 )/ 2 | J B'sin X '|^+ 1 )/ 2 sin2x, 
V = 0, (B7) 
where we reinstated the cosmological factor 1 + z. The function is 



V3 r /^3p-l\ ^3p + 7 



e 



3 



3e 



27rm^c 5 j 



(p-i)/2 



I/ -(p-i)/2 2fej (B8) 



4 V 12 / V 12 / m eC 2 

where e and m e are the charge and mass of an electron, is the Euler gamma-function. 

The degree of polarization of the observed radiation pulse is expressed as II = sj Q 2 + U 2 /I, 
giving Eq. (2). The resultant position angle of the electric field Xres measured by the observer is 
found from 

cos 2 * rcs = 7^2=p ' sin 2 * rcs = 75i=p ' ^ ^ <7r - (B9) 

It can be checked that in our case under the change of <p to 2ir — <p in the integrals (B7) the value 
of Q is not changed, and the sign of U is reversed. Therefore, the Stokes parameter U integrates 
out to zero. Consequently, if Q > then x rcs = 0, if Q < then Xms = Thus, the observed 
electric vector can be either parallel or perpendicular to the projection of the axis of the flow on 
the plane of the sky. For a shell carrying only toroidal magnetic field Q > 0. 



C. Lorentz transformations of polarization vector 

In this appendix we first derive Lorentz transformations of polarization vector of the linearly 
polarized radiation emitted by a relativistically moving plasma with a given magnetic field and 
then find an angle x between a given direction (chosen later as a direction along the projection of 
the axis of the flow in the plane of the sky) and the direction of linear polarization of the waves for 
a spherically expanding shell. 

Let n' be a unit vector in the direction of a wave vector in the plasma rest frame, B' be a unit 
vector along the magnetic field in the plasma rest frame. The electric field of a linearly polarized 
electromagnetic wave is directed along the unit vector e' = n' x B' and the magnetic field of the 
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wave is along the unit vector b' = n' x e', such that the Pointing flux along e' x b' is directed along 
n'. We will make a Lorentz boost to the explosion frame to find an electric field e there, normalize 
it to unity, and project e on some given direction (e.g. , along the projection of the flow axis on 
the plane of the sky) . 

Fields in the wave expressed in terms of the direction of a photon in the explosion frame n are 

„ n x B' l + r(l-n-v) 

p! _ i v —TV y v 

r(l-n-v) (l + r)(l-n-v) ' 
CJ ty , ( B ' " l+r(l-n-v) p \ 

^ ^r^i-n-v) 2 r(i + r)(i-n-v) 2 J 

I / (i + r(i-n-v)) 2 p l + r(l-n-v) p \ 

+ ( v( i + r)2(i- n .v)2 B ■ v -r(i + r)(i-n.v)2 B ' n ) v - (C1) 

It may be verified that e' x b' is still directed along n'. 

Next, make a Lorentz transformation of e' back to the lab frame 

e = r(e'-— !— (e'-v)v-vxb') (C2) 



r + i 

and normalize to unity. After some rearrangement we find 



n x q' 



vV-tn-q') 2 ' 

q' = B' + nx (vxB')-^p(B'-v)v. (C3) 

Finally, we may express the rest frame unit vector B' in terms of the laboratory frame unit 
vector B. Assuming ideal MHD, there is no electric field in the rest frame of the plasma, E' = 0. 
Then, we obtain 

B = , 1 (& - ^(B' • v)v 



1 - (B' • v) 2 

i + r)B + 
(i + r) v / i + r 2 (B-v) 2 



B' = (l + nB + r 2 (B.v)v 



to get 



n x q 

e = 



y/q 2 - (n- q) 2 

q = B + nx(vxB). (C5) 

This is a general expression giving the polarization vector in terms of the observed quantities B, n 
and v. If, for a moment, we adopt a frame aligned with the direction of motion (Fig. 4), we find 
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from (C5) 



cos(0 + ib) — 8 cos ib _ N 

tang = cot7/ ~ ~ ^, C6 

1 — p cos 9 



reproducing Eq. (16) in Blandford & Konigl (1979). 

In our particular case B' • v = 0, so that the fields in the rest frame of the emitting plasma 
element and the laboratory frame are aligned B' = B. The general relations then simplify. Setting 
B' • v = in Eq. (C3) gives 

n x (b + (n x (v x B))) 

e = . V V 11 . (C7) 



- n • v) 2 - (B • n) 2 /r 2 



Next we introduce a unit vector 1 normal to the plane containing n and some given direction 
(in our case the direction of the projection of the axis of the flow to the plane of the sky). Then, 

cosx = e • (n x 1), sinx = e • 1. (C8) 

Using Eq. (C7) we find 

cos x = 1 ( B • ( 1 + ^ ' V ^ n 



^-(nfei) ; 



1 — n • v 



r /. A , r(n-v)(B-vxl) T(B-nx v)(l-v) 
sin x = —— — I (B • 1 x n) H A '- + 



1 



B n 



r(l-n-v) 



2 \ v ' (i + r) (i + r) 



In our case 

B = {— sin (p, cos (b, 0}, n = {sin# b, 0, cos 6> t>}, 
v = /3{sin# cos 0, sin sin<^>, cos 0}, 1 = {0, 1, 0}, 

n • v = P/j,, jj, = cos G = cos cos 6> Q b + sin 6 sin o b cos 4>, (C10) 

which gives Eq. (7). 



D. The ultra-relativistic limit. 

In the ultra-relativistic limit, when r » 1, the maximum polarization fraction for a given 
velocity field may be found analytically. Because of the Doppler boosting effect described by the 
factor T>, the contribution to the integrals in formula (2) comes from the small patch <f> ~ 1/T and 
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Q — #ob ~ 1/r. We can introduce rescaled variables £ = Tip and ^ = T(9 — ob ) and change the 
integration from 9 and <p to £ and ^. The integration limits can be taken from — oo to +00 for 
both £ and Making expansions for T S> 1, £ ~ 0(1) and * ~ 0(1) the expression (7) results in 

cos 2v = (^-^^ 2 ^ + l) 2 -4^ 2 sin 2 g ob 

COSZX (*2_ C 2 sin 2^ ob+ l)2 +4 x|>2£2 sirl 2 ob > ^ 



the expression for sin x' becomes 

,2 _ 1 

(l + * 2 +£ 2 sin 2 # ob ) 2 ' 



• 2 / n 4£ 2 sin 2 fl ob 



and the expression for T> becomes 



IT 

D = — (T)3) 

^2 + £2 sin 2 + I ■ V ! 



Note, that the variable £ and sin# ob enter in the integrals (2) only in the combination £sin# ob . 
Therefore, by changing integration to a new variable £1 = £ sin # ob the value of the integrals and IT 
becomes independent of # ob . Therefore, for T S> 1 the polarization degree is insensitive to # ob as 
long as # ob » 1/r (see Fig. 3). 

Further, it is convenient to switch to the integration in "polar" coordinates a and r in the plane 
£1 and ^, which are introduced according to a = £ 2 + \I> 2 , £1 = \fo cos(r/2), and ^ = v / o : sin(r/2). 
After some algebra we obtain 



roc r2n 

/ dcr / dr(l + <r 2 cos 2r — 2cr cost) 
■/ ■/ 



(l + a 2 -2CTcosr)(f- 3 )/ 4 



n = Jt±Lk k : : (i±£g ( D4) 

P + 7/3 /■«> f 2 - (l + a 2 -2acosr)(f+ 1 )/ 4 k j 

/ a<7 / GST rrr- 

Jo Jo (l+a) 2 +P 

For p = 3 expression (D4) gives IT = 9/16 « 56%. This value is the value of the horizontal 
asymptotic of r(# ob ) curve for p = 3 in Fig. 3. For arbitrary p asymptotic values of the polarization 
are plotted in Fig. 5. 
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Fig. 1. — Geometry of the model. A narrow shell Ar ~ t s c ~ 3 x 10 cm, dominated by the 
toroidal magnetic field, expands quasi-spherically with angle dependent velocity v(0). The observer 
is located at an angle 6> b with respect to the polar axis. 
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Fig. 2. — Polarization map of the emission from GRB on the observer plane of the sky in rectangular 
coordinates I and s. Axis s is directed parallel to the projection of the axis of the shell on the sky. 
The observer line of sight makes 30° angle with the axis of the shell, # b = 30°. Plots are made for 
four different values of T and p = 3. Solid circles have radii 1/r. As the intensity of radiation is 
highly peaked in the area of the size ~ 1/r, we zoomed in on this area in the plots with f>l. 
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Fig. 3. — Dependence of the polarization fraction II = Q/I on the viewing angle # b f° r different 
Lorentz factors T for isotropic expansion (T(6) = constant = 10,50,100, right to left) and the 
power-law particle distribution dn/de = e~ p ; upper curves p = 3, lower curves p = 1 (for p = 2 the 
asymptotic value is 43%). At o b = polarization is zero, growing to large values when 9 Q b > 1/r. 
Depolarization of emission due to differential rotation of the position angle of the linear polarization 
in the n — v plane reduces the maximum possible polarization fraction below the theoretical limit 
for a homogeneous field. 
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Fig. 4. — Characteristic swing of polarization angle due to relativistic motion. In the frame aligned 
with v, the electric field of the wave e and the observed magnetic field B make angles £ and ij 
with the plane containing v and n, while their projections make angles G and ip with the axis z 
perpendicular to v — n plane (after Blandford & Konigl (1979)). 
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Fig. 5. — Asymptotic value of polarization for r — > oo as a function of p for spherically divergent 
flow (solid line), maximum polarization for a homogeneous stationary magnetic field Il max (dashed 
line) . 



